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I. INTRODUCTION 



Many efforts in the literature, over several decades, have been devoted to the construction of gauge theories of 
gravity, mainly motivated [l[ by: 



o 

I 

5h 

(i) conceptual simplicity within the geometric framework of fibre-bundle picture of modern field theories; 

I ■ (ii) the search for quantum gravity (see Q and the references therein) , that should unify quantum theory and general 
\ relativity. Such an unification is logically compelling if one thinks that gravity couples to the energy-momentum 
tensor of matter, but matter fields are defined and treated via quantum theory, and the semiclassical analysis where 
' gravity remains classical is a hybrid scheme that cannot hold at all energy scales. 
' ... 

• (iii) attempts of building unified theories of all fundamental interactions. The various models can be classified 
' according to the relevant structure group in the fibre jl|. Within the Poincare group, this has been either the 
translations, or the Lorentz group, or the whole Poincare group. The enlargement of the group has involved using the 
conformal, or special affine or affine groups. Moreover, the so-called supergroups have included the graded-Poincare 
or the superconformal group. Interestingly, re-expressing the action of supergravity in Yang-Mills form has made it 
possible to simplify enormously the proof of its invariance under local supersymmetry transformations [|[ . 

An alternative view is also available, which is convenient, if not natural, if one is interested in the geometry of the 
space of histories and its application to functional integrals in quantum field theory Q ■ Within this framework, gauge 
theories are all field theories whose action functional S is invariant under an infinite-dimensional Lie group. Gauge 
invariance is then expressed by the existence of vector fields Q a (a being a Lie-algebra index) on the space of histories 
that leave invariant S, i.e. 

Q Q S = 0. (1.1) 

All known gauge theories belong to one of the following families: 

(i) Type I: The Lie bracket [ , ] of such vector fields is a linear combination of the Q a themselves, i.e. 

[Qa,Qp] = Q-y, Cip = -Cj a , (1.2) 
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the coefficients CLg being structure constants, hence independent of the field variables ip 1 . 

Interestingly, Maxwell's electrodynamics, Yang-Mills theory and Einstein's gravity share the property of being type- 
I theories, which was overlooked for a long time and even nowadays, since Yang-Mills theory was first formulated in 
Minkowski space-time, and only there can its perturbative renormalizability be proved, this property being spoiled 
by space-time curvature [J]. 

Type II: The structure constants of the previous equation are replaced by structure functions C 1 a * , i.e. dependent on 
field variables, so that 

W dl ^°- (L3) 
An example is given by simple supergravity in 4 space-time dimensions, endowed with auxiliary fields. 

Type III: The right-hand side of the defining equation acquires also an inhomogeneous term, so that 

5S 

1 Sip' 



[Q a ,Qf>] = cr a ^ + Ti —, (i.4) 



where T is skew-symmetric in the Lie-algebra indices. An example is provided by simple and extended supergravity 
in any number of space-time dimensions, without auxiliary fields. 

Now we can revert to a topic more immediately related to the calculations we are interested in, i.e. the formulation 
of Einstein's theory as SL(2,C) gauge theory, following We study in detail the equations for nonlinear gravito- 
electromagnetic waves and their solutions. 



II. ELECTRODYNAMICS AND GENERAL NONLINEAR CONSTITUTIVE EQUATIONS 

To set the context and specify notation, we review established results in this section We consider only flat 
spacetime, so that the metric is given by the Minkowski tensor rj^v = (1, —1, —1, —1), x p = (ct, x l ), fi, u, . . . = 0, 1, 2, 3; 

— 1,2,3; with d p = d/dx p = [c _1 c)/c?i, V] . The antisymmetric Levi-Civita tensor is written e pvpa ', with 
£ oi23 _ with tensor notation, we denote by the electromagnetic field tensor and by G M „ the tensor expressing 
the constitutive equations (see below), while their Hodge dual tensors are F pv = ^e puf " y F pa and G pv = pa G pa , 
and hence Maxwell's equations become 

d^ v = 0, d^ u = f, (2.1) 

where = (cp, j) is the 4-current. The first equations in (|2.1[) imply F pv = d^A v — d v A„, where A p is an Abelian 
gauge field; but in general there is no similar representation for G pv . With covariant notation we have 2 independent 
invariants 



A = -F^F»\ Y = \F^ V F^. ( 2.2 ) 



Then the nonlinear constitutive equation of @, H| is 

Gfj, v = N (A, Y) F M „ + cM (A, Y) F^ . (2.3) 

Although the constitutive equations (I2.3[) are fairly general, they do not take account of a variety of possibilities, such 
as anisotropic media [9(, chiral materials where derivative terms enter [10(, piroelectric and ferromagnetic materials, 
and so forth. Therefore it was proposed Q to generalize (|2.3I) by introducing the constitutive tensors S^, R p plJ , and 
1 j ft = 1, 2, 3, . . . and write a general nonlinear constitutive equation 

Gfi, v = V + RP° v F pa + QP^D Xl F p<x (2.4) 
+ QP^D^D^Fpv + ... + Q^ - x "D Xl D X2 ...D Xn F pa . (2.5) 

Evidently the formula (|2.4p . taken together with Maxwell's equations, includes all possibilities discussed so far, 
as well as new ones, defining the general nonlinear electromagnetic theory. We impose Lorentz covariance on the 
constitutive equations, the constitutive tensors will depend on the fields through the invariants A and Y", as follows: 
Sft.,/ is a constant independent of A and Y , while 

R%=R%,(X,Y), Qfru Xl "' Xn =Q^ Xl - Xn (X,Y,...), (2.6) 
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where ". . . " denotes covariant derivatives of the invariants X, Y up to n-th order. Obviously Sp V is antisymmetric, 
R p ° is antisymmetric in its upper and lower indices separately, and the Q p ^ v x " n are antisymmetric in their upper 
and first two lower indices; with respect to the Ai, they are symmetric Lorentz tensors. 

If the equations of motion for the nonlinear theory derive from a Lagrangian L (X, Y) which is a scalar function of 
the invariants X and Y but does not depend on their derivatives, then from the usual definitions together with Eq. 
p.3p we have 

_ dL(X,Y) dL(X,Y) ~ 

Comparing (|2.4[) and (|2.7I) gives us the constitutive tensors, 

S^u = 0, 

dL(X,Y) dL(X,Y) \p 5a 

^ ~ dX ^ v] dY £ ^ xsr] 11 ' ( 2 - H > 



Q p.a 1 ...a„ = 

Now we can relate this formulation of nonlinear electrodynamics with gravity. 



III. GRAVITY AS A GAUGE THEORY 

Here we consider the formulation of the Einstein gravity as SL (2, C) gauge theory [f|. The gauge potentials are 
interpreted in terms of an afiine connection on a complex vector bundle. In the SL (2, C) gauge theory they are taken 
to be the dyad components of the spinor connection, and are therefore the Newman-Penrose spin coefficients. 

Let r B be the complex components of the spinor connection V, i.e. 

v^/ = r/ B e a s , (3.1) 

where (™ is a spinor basis of the complex vector bundle. The dyad components of the spinor connection are defined 
by 

B p a b = r M A B £a A t b B- ( 3 - 2 ) 

The field strength of the gauge theory (spinor curvature tensor) can be expressed through r AB as 

Fpv A B — dyT^ A B +^ t iA C ^vC B ~~^vA C ^ l iC B - (3-3) 

Thus, we can rewrite F A B directly through the Riemann curvature tensor R a ji^ v as follows: 



where er Qa &' = <T aAB'£, c ^£,b' > the <r a AB' being the Infeld-van der Waerden symbols which turn tensors into spinors 
and express the isomorphism between the tangent space at a point of space-time and the tensor product of unprimed 
spin-space and primed spin-space. In Minkowski space-time they reduce to the Pauli spin matrices. 
The field strength obeys the Bianchi identity 

V P F a B + permutation = 0, (3.5) 

which are indeed the second set of the Maxwell equations. The Riemann tensor R a pp.v can be decomposed into the 
Weyl tensor W a pp V , the Ricci tensor R^ v and the Ricci scalar R as 

Rp,vpa — Wfiupa + gp,[pR a ]v + 9v\pRa\p, + ^Rg fi[<j 9 p\v ■ (3.6) 

The Einstein equations are 

R^v = k \ Tp lV - \gp V T ) . (3.7) 
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Thus, the decomposition (|3.6[) becomes 



Rfiupa — W pupa + 2/c5[ p [ M i?i / ] (T ]] - -kg p [ a g p ] u T. (3.8) 
A corresponding expression can be written for the field-strength tensor 

F^v = F^ v + kr^u, (3-9) 

Ffiv = G^ v — Sfj, u (3.10) 

Gfj, u = Sfj, u + R p ° v F pa + =$> (3-11) 

Einstein gravity R% = 5ffi, S pu = -fcr^, Q = 0. (3.12) 

This suggests the general structure 

S llv =C{x)5l5l ] e [llv]X5 , (3.13) 

i?£ = A (X, Y) 6^ + B (X, Y) e lfiu]xs v Xp V 6a , (3-14) 

while Q is an as yet unspecified tensor Q P °J^ ■ 
Let us introduce the covariant derivative 

D^X = d fl X-[B^,X], (3.15) 

then we can write the Bianchi identity as 

e w D {B)P pa = o. ( 316 ) 

By using the formula (|3 .0[) we obtain 

£ ^P° D {B)PW _ fcJM ( 317 ) 

where J M is the current corresponding to the matter 

Now we are ready to compare the above formulas with the standard nonlinear electrodynamics: the Maxwell 
equations (|2.1[) and the constitutive equations (|2. 41) . We can identify (|2.4|) with (I3.9|) and obtain 





- pw 

x (11/ ' 


(3.19) 


F 


— p 


(3.20) 




kT pv , 


(3.21) 


f 


= kj p . 


(3.22) 



It follows from (|3~TTj) that 

£ w D {B) Gpa = f. (3.23) 

To obtain the analogy with the first set of the standard Maxwell equations d^G^ = j u , we introduce the Hodge 
dual tensor 



1 

2 

Then we have 



G pv = -e^ pa G pa . (3.24) 



D (B)^u = }_ f (3 25) 
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IV. NONLINEAR GRAVITO-ELECTROMAGNETISM 



Here we will turn to 3-dimensional notation for convenience and examples. Let us consider the gravito- 
electromagnetism in the weak- field approximation (following, e.g., [HI). Recall the standard Maxwell equations 
in SI units [H[ 

dB 

curlE = — — , divB = 0, 
at 

<9D 

curlH = — +j, divD^p, (4.1) 
at 

where E is the electric field, B is the magnetic field, p is charge density, j is electric current density. In the linear case 

B = p H, D = e E, (4.2) 
In the nonlinear case these equations can be presented in the form [l3T ] 

D = M (h,I 2 )B + \n (h , I 2 ) E, 

H = JV(7i,/ 2 )B-M(7i,/ 2 )E, (4.3) 

where the invariants are 

I 1= B 2 -\e 2 , I 2 = BE. (4.4) 

&■ 

f)B 

curlE g = — divB ff = 0, (4.5) 

1 1 1 

curlB 9 = ? -^ + — divE g = -p 9 , (4.6) 

where E g is the static gravitational field (conventional gravity, also called gravitoelectric for the sake of analogy), B g 
is the gravitomagnetic field, p g is mass density, j ff is mass current density, G is the gravitational constant, e g is the 
gravity permittivity (analog of £o)- Here 

M s = 2~, (4-7) 



Their gravitational analogues in SI are 



In the linear-gravity case 



are the gravitational permittivity and permeability, respectively. 

The main idea is to introduce analogues of H and D to write (|4.6I) in the Maxwell form for 4 fields in SI as 

f)B 

curlE 9 = -^, divB s = 0, (4.8) 
curlH g = ^+j s , divD ff = Pg . (4.9) 

D 9 = e 3 E s , (4.10) 
B 9 = p g H gi (4.11) 

£ g »g = ^- (4-12) 

The main idea is the following: the linear-gravity case (|4.10[) - (14.12[) corresponds to weak approximation and some spe- 
cial case of gravitational field configuration. We generalize it to nonlinear case which can explain other configurations 
and non-weak fields, as in (|4.3p . by 

T> g = Mg (l g l,I g2 ) Bg + ^Ng (Ig^Itf) Eg, (4.13) 

H 9 = N g (I gl ,I g2 )B g - M g (/ Sl ,/ 92 )E 9) (4.14) 

where the invariants are 

hi = V 2 g - 4^, I g2 = B g ■ E g . (4.15) 

The gravity-Maxwell equations (|4.8p ~ (|4.9p together with the nonlinear gravity-constitutive equations (|4.13[) - (|4.14[) 
can give a nonlinear electrodynamics formulation of gravity (at least its particular cases). 
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V. LINEAR GRAVITO-ELECTROMAGNETIC WAVES 



curlE 9 = — 5 /,divB s = 0, (5.1) 
curlH g = ^,divD g =0 (5.2) 



The gravity-Maxwell equations for gravito-electromagnetic waves (far from sources) are 

dB 

dT> 

with generic values of permittivity and permeability (14.71) . Then 

curl E 9 = - Mg ^,divH g = 0, (5.3) 
curlH g = £g ^,divE g = 0. (5.4) 

d <9 2 H 1 9 2 H 

We differentiate the first equation with respect to time: curl — E g = — fi g 9 => — curl(curlH g ) = — /x g 9 . 

Since curl (curl H g ) = grad(divH 9 ) — AH 9 = — AH g , then 

~~di 

n _ <9 2 E„ -1 , , s d 2 E 

9 ~ £ 9 

equation for E g , 



AH S =£ A ^/. (5.5) 
By analogy, from the second equation curl T-H g = e g " 9 =4> — curl (curlE g ) = e g 9 . Hence we get the wave 

C/V CJX/ f^Q C/~C 



VI. NONLINEAR GRAVITO-ELECTROMAGNETIC WAVES 

The differences begin with the constitutive equations f!4. 13)) (|4. 14[) . For simplicity put M g = 0. Then 

D s = ^E g , (6.1) 



1 

N 

where N = N g (7 gl , I g 2)- The Maxwell equations become (hereafter the dots denote time derivatives) 

IV H -1^ 

n 9 n at 



B s = ^H g (6.2) 



curlE S = -U? H 9 --^> (6-3) 



(r J \c 
Take derivative of (|6.3[) with respect to time and get 



' »v ( ^H g j - H g grad (I J + I div (H g ) = 0, (6.4) 
N N <9E 

curlH g = -E g + -^, (6.5) 

div ( -E g = E g grad - + - div (E 9 ) - 0. (6.6) 



(9E c 2 N 

From (16.51) it follows n 9 = — curlH„ E„. Then we get 

v dt N 9 N 9 



curl ^curlH g - I Eg = - ( I )" H g - 2 ( I V & - 1^. (6.8) 
I AT 9 iV 9 j \ AT/ 5 \N dt N dt 2 v; 
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The left-hand side here is 

(c 2 N \ 

curll-curlH ff --E fl J 

c 2 c 2 c 2 N N 

= grad — x curlH s + — grad div H g - -AH 9 - — curlE g - grad — x E 9 . 

From (16. 4p we get div (H s ) = — JVH g grad (jzj 0- Thus, the nonlinear analogue of the wave equation is 



c 2 c 2 c 2 N N 

grad — x curlH 9 + — grad div H g - -AH 9 - — curlE s - grad — x E g (6.9) 

1 \ " ( 1 \ ' <9H„ 1 d 2 H a 

H 9 - 2 7F -7T V-T77 2 - (6-10) 



4 JVy tf \N J dt N dt 2 

Note that if iV = const, then we obtain the usual wave equation 

1 <9 2 H 
^~dt 2 

Take now the constitutive equations in the form 



AH 9 = -^/. (6.11) 



Ig — ±1 J_.g 

where JV, M are constants. In absence of sources, the Maxwell equations become 

dB 



N 

D 9 = MB g + ^E s , (6.12) 

H 9 = NB g - MEg, (6.13) 



curlE s = -^£, divB g = 0, (6.14) 



<9D 

~di 

If we express the Maxwell equations through E 5 and B g , then second pair of equations become 

"a* 



curlH g = — ^, divD g = 0. (6.15) 



curlHg = ^/, (6.16) 



JV curlB, - M curlE, = M— ■£ + (6.17) 

at c 2 dt 



dB 

Since curl E ff = — , from the last equation we get 

1 <9E 



curlB s = -— f. (6.18) 
d- at 

N 

The second equation, div D g = 0, reduces to M div B ff + — div E g = 0. Since div B 5 = 0, we get 

divEg^O. (6.19) 
Thus, using constitutive equations with constant M and JV we have Maxwell equations in terms of B g and E fl , i.e. 

OB, 

dt 



curl Eg = divBg = 0, (6.20) 



1 <9E 

curlBg = -2^, divE g = 0. (6.21) 
At this stage, we get the wave equations in the standard way. The time derivative of the first equation yields 



<9„ <9 2 B„ „ ,„ s 5 2 B 

-. Since 

1 <9 2 B 



curl^-Eg = - - ' y => c 2 curl (curl Bg) = — ~q^T- Since curl(curlBg) = grad(divB g ) — AB 9 = -AB S , then 



ABg = ^. (6.22) 
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d 1 c^E 1 <9^E 

By analogy curl — B, = - „ g =>■ — curl (curlE„) = — 9 , and we get the wave equation for E„, 
at cr ct cr ot 2 

AE S . 1^. (6.23) 

Thus, the gravi-electromagnetic waves E g and B g have speed c and do not depend on the constants M and N. 

VII. WAVES AND CONSTITUTIVE EQUATIONS FOR LINEAR CONSTITUTIVE FUNCTIONS 

Let us consider the constitutive equations (|4.13l) - (|4.14l) as linear functions of invariants, i.e. 

M = M g {Igulgi) = a m I g i + b m I g2 , (7.1) 
N = N g (I g i,I g2 ) = c 2 e g + a n I gl + b n I g2 , (7.2) 

dm i b m ,a n , b n being some constants. From all the Maxwell equations in material media, and in the absence of sources 

one finds curlH„ = — ^ curl(7VB„ - ME„) = — MB„ + — E„ ) , and TV curl B„ — M curl E„ = M— ^ + — — ^. 

dt a dt \ c 2 J dt c 2 dt 

<9B 

Since curl E g = — , from the last equation one gets 

crlB^i^. (7.3) 

( N \ N 

The second equation, div D g = 0, reduces to div I MB 9 H — jE s I = 0, or M div B g + — div E g = 0. Since div B g = 0, 

one gets 

divE g = 0. (7.4) 

VIII. INVERSE PROBLEM OF NONLINEAR GRAVITO-ELECTROMAGNETISM 

In electrodynamics the direct solution of the Maxwell equations together with the nonlinear constitutive equations 
is a nontrivial and complicated task even for simple systems 0, H|. In previous sections we presented some very 
special cases of the nonlinear functions ./V and M. Here we formulate the following inverse problem: if we have some 
particular solution of the gravity-Maxwell equations (|4.8[) - (l4.9p . can we then find the exact form of the corresponding 
nonlinear gravity-constitutive equations (|4.13[) - (|4.14[1 ? 

It is natural to consider the case of plane gravitational waves, when the fields have only one space coordinate. We 
will show that even in this case one can have a nontrivial nonlinearity. Let us choose E g and B g mutually orthogonal 
and perpendicular to the direction of motion 





E 9 = y j , Bj=l II I . (S.I) 

where E = E (t,y), B = B (t,y). Now the invariants (|4.15[) become 

Igi = B 2 - \e 2 = I, (8.2) 

cr 

Ig2 = 0. (8.3) 

The use of the nonlinear gravity-constitutive equations (|4.13[) - (|4.14[) gives for other fields 

/ 1 



D 9 



\ MB 




NE 

I » H fl = | 1 . (8.4) 
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where N = N (I), M = M (I) are the sought for gravity-constitutive functions. They depend on 7 only, because 
of Lorcntz invariance (see @, H|). Inserting the fields (|8.1[) and (|8.4[) into the gravity-Maxwell equations (|4.8[) — (|4.9[) 
without sources gives us 3 equations (hereafter, a prime with the corresponding subscript denotes the first partial 
derivative with respect to the variable in the subscript, while dot denotes time derivative) 

E' y = B, (8.5) 

(NB)' y = ^{NE)\ (8.6) 

(ME)' y = (MBy . (8.7) 

Now we take into account that the gravity-constitutive functions N, M depend only on the invariant / and present 
T6|) - (|8.7I) as the differential equations for them 

N 'i ( BI v - ^ Ei ) + N ( B v- ^ E ) = °> ( 8 - 8 ) 

Mr 



where we have exploited the identities 



r; (EI' V - Bl) = 0, (8.9) 
N' y = N'jl'y, M' y = M'jl' y , (8.10) 



N = N'jl, M = Mji. (8.11) 
The second equation (|8.9|) can be immediately solved by 



M{I) = { A/ °= const < sEr v *Bi, 

1 arbitrary, if EI' y = BI. K ' 



The first equation (18. 8[) can be solved if 



( B v f) 

A = 7^— 4- (8.13) 

{bi> - f ) 

depends only on /, which is a very special case. One then has the differential equation 

TV; + A (!) TV = 0, ( 8 _ 14 ) 

and its solution is 

N(I) =^ c" /A(/)d/ . (8.15) 
Otherwise, by using the expressions for I' and I from (|8.2p . i.e. 

4 - 2BB v -^r"'^ 21313 -—' ( 8 - 16 ) 

we obtain 

2N 'i ( B2B v + ^ E ' E ^ EBE 'v) + N ( B 'v- i> E ) = °» ( 8 - 17 ) 

where the sum of terms in brackets is not a function of I in general. 

Usually, in the wave solutions the dependence of fields on frequency w and wave number k is the same, and therefore 
we can consider the concrete choice 

E(t,y) = f(eut + ky) = f(X(t,y)), B (t, y) — g (suit + ky) = g(X(t, y)), (8.18) 
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where s = ±1, with / and g arbitrary smooth nonvanishing functions. Bearing in mind that 

E' y = fxK = kf x , B' y = g' x X' y = kg' x , 



our Eq. (|53)1 yields 



Therefore 



E = f x X = euf' x , B = g' x X = eug' x , 

k .f'x = ^g'x- 



>(X) = —f(X) + a, 

SUJ 



(8.19) 



(8.20) 



where a is a constant, so that both E and B can be expressed through one function only, i.e. /, and the invariant I 
reads eventually as 



I=\[ k 2 - =g )f + 2— af + a 2 

UJ Z \ C z J EU) 



The equations for the gravity-constitutive functions take therefore the form 



Ni 



21 k 2 



LU 2 \ U 2 , 



+ N[k 2 ~— ) =0, 



Mil 



X 



V I U2 



SLC 



2ka 



a = 0, 



where we have exploited the identities 



p fl / ,/ // 



suj 



3.21) 

3.22) 
3.23) 

3.24) 



fjf, 



I f f fj\ 



X 



./' ( /.-' ) k:^c<\ 



and, after some cancellations, 



/ ( k j ) + fcewa 



2 fc 



2/ 



fc 2 - 



2 — a 

SUJ 



UJ~ 



1 + 2 — a', 



while 



f% - gi = iff' - gf) 



2f 



n k 
2 — a 

SUJ 



(8.25) 



(8.26) 



.27) 



ffL - gf = .f'x( k f - £uj g) = -euf x a 



(8.28) 



The results of our analysis now depend on whether or not a vanishes. Indeed, if a = 0, M is arbitrary and hence 
we obtain the equation 



(2IN'j + N) = 0, 



3.29) 



which implies that either the dispersion relation 



fc" - = 



(8.30) 
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holds, with N kept arbitrary, or such a dispersion relation is not fulfilled, while N is found from the differential 
equation 

2IN'j + N = 0, (8.31) 

which is solved by 

N(I) = ^. (8.32) 

By contrast, if a does not vanish, M equals a constant Mq, while N solves the more complicated equation (|8.22l) . 
At this stage, to be consistent with the dependence of N on / only, we have to require again that the dispersion 
relation (|8.30[) should hold, jointly with N[ = 0, which implies the constancy of N: N — Nq. 



IX. CONCLUDING REMARKS 



We have proposed a way to take an extremely general approach to a nonlinear formulation of gravity as classical 
electrodynamics. The framework formally includes nonLagrangian as well as Lagrangian theories, and accommodates 
the description of nonlocal effects. This is accomplished through generalized constitutive equations (and constitutive 
tensors). We expect future research directions to include the detailed development of new examples within this 
framework. 
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